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We  develop  a  modified  Ox 8  envelope-function  approximation  (EFA)  formalism  for  the  noncommon-atom 
(NCA)  superlattices  (5L  s),  incorporating  the  effect  of  anisotropic  and  other  interface  (IF)  interactions  that  go 
beyond  the  standard  EFA,  The  boundary  condition?  in  Lhe  presence  of  IF  interactions  ufe  used  to  set  up  a 
secular  equation  (including  r  transfer  matrix  derivation)  whose  physical  transparency  makes  possible  a  number 
of  valuable  insights  (possibility  of  IF  bound  states,  analytic  solutions,  Indirect  gaps,  elo.).  We  show  that  the 
heavy -holc-spin-orth  IF  coupling  is  very  important  due  to  die  IF  localization  of  the  SO  wave  function 
components  and  the  ability  of  the  IF  potential  to  potentially  bind  a  hole  at  the  IF's,  all  of  which  pose 
convergence  problems  for  perturbative  solutions.  With  two  adjustable  parameter  for  the  two  possible  IF’s,  we 
Find  a  very  good  agreement  between  experiment  and  theory  for  the  band  gaps  of  several  sets  of  vety  Song- 
infrared  and  midinfrarod  InAs/GaSb  SL’s  grown  at  sevcrnl  laboratories  and  by  us.  The  band  gaps  as  a  function 
of  GaSb  and  ItiAs  widths  axe  explained  in  terms  of  variations  of  the  l  IH  and  condition  (C)  band  bandwidth*. 

We  show  (hat  the  cut-off  wavelengths  can  be  reduced  by  increasing  the  GaSb  layer  width.  Thus,  a  consistent 
application  of  the  EFA  method  with  the  inclusion  of  well  established  IF  effects  can  provide  useful  physical 
insights  and  possesses  good  predictive  capacity  in  the  denigu  of  NCA  SL’s. 
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L  INTRODUCTION 

InAs/InGaSb  superiattices  (SL's)  have  been  proposed  by 
Smith  and  Mailhiot!|3  and  grown  by  many  groups  for  use  as 
detectors  and  lasers  operating  in  the  short  to  veiy  long  infra¬ 
red  wavelength  regions. 5-5  This  5L  system  has  been  modeled 
using  a  variety  of  computational  approaches,®-33  including 
the  popular  envelope-function  approximation  (EFA)  method. 
Recently,  however,  problems  wilh  she  numerical  accuracy  of 
the  standard  EFA  in  the  case  of  InAs/GaSb  binary-binary 
superlattices  have  been  reported.*'1,1'21  For  example,  North¬ 
western  Sample  No.  1303 — (InAs)Lj/(GaSb)4  with  InSb 
IF's  (Ref  12) — was  measured  to  have  the  cutoff  of  about  30 
whereas  our  standard  8X8  EFA  code32  predicts  a  5.5 
pirn  cutoff.  Another  Northwestern  SL — (LnAs)|7/(GaSb)T 
with  InSb  IF's— has  the  cutoff  of  19  jum,n  whereas  our  stan¬ 
dard  8X8  EFA  predicts  a  7.5  ^tm  cutoff.  In  our  EFA  simu¬ 
lations,  changing  the  layer  widths  by  one  monolayer,  varying 
band  offsets  by  ±10  meV,  and  even  relaxing  the  SL  did  not 
change  the  cutoff  by  more  than  1  jrm. 

The  standard  EFA  encounters  the  same  problem  for 
mrd-IR  IriAs/GaSb  Si's.  For  example,  Kaspi  et  al 23  grew  a 
series  of  SLTs  and  measured  their  absorption  cutoffs  as  a 
fimetion  of  GaSb  width.  In  disagreement  with  Kaspi ns  data, 
the  standard  8X  8  EFA  model  does  not  predict  large  changes 
in  the  band  gap.  Whereas  the  EFA  calculated  band  gaps  for 
InAs/GaSb  SL+s  varying  from  6  ML/6  ML  to  6  ML/24  ML 
increase  from  421  to  452  meV,  absorption  and  PL  data  indi¬ 
cate  that  the  gaps  increase  more  rapidly  from  about  320  to 
430  meV22,23  The  larger  error  for  the  thinner  SL  suggests  a 
correspondingly  greater  role  of  the  IF's.  Similarly,  in  the  se¬ 
ries  of  InAs/GaSb  SL  samples  ranging  from  8  ML/3  ML  to  8 
M1740  ML,  the  band  gaps  via  the  standard  EFA  range  from 
343  to  352  meV,  whereas  absorption/PL  data  range  from 
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(277/263)  to  (352/353)  meV,22  with  disagreement  again  in¬ 
creasing  for  thinner  SL’s. 

This  disagreement  motivated  us  to  reexamine  the  standard 
EFA  model. J*34  Recently,  the  standard  EFA  has  been  ex¬ 
panded  to  model  IF  anisotropy  in  NCA  SL's  (Refs.  35-39) 
and  used  to  account  for  (heir  giant  absorption  anisotropy.  In 
HI  V  heteroj unctions  with  noncommon  anions  or  cations* 
grown  in  the  [00 i]  direction,  adjacent  rows  of  atoms  are  at 
right  angles  to  each  another,  so  that  the  true  IF  symmetry  is 
CZl,  rather  than  £>2J  as  assumed  in  the  standard  EFA.3*  The 
lowering  of  the  symmetry'  causes  heavy-hole- light-hole 
(liH-LH)  interaction  even  at  the  center  of  the  Brillouin  zone 

(BZ),  point  T.  Using  the  theory  of  invariants,40,41  I  vchenko, 
Kaminski,  and  Rossler42  showed  shat  she  form  of  the  IF  in¬ 
teraction  is  the  consequence  of  die  symmetry  of  SLFs  with 
NCA  IF's  and  modeled  it  with  a  short-range  delta  funelion 
potential  centered  at  the  IF's.42-4*  Foreman47  than  formally 
justified  the  use  of  short  range  IF  potentials  to  model  IF 
effects  at  NCA  IF’s  using  Burths  EFA  representation  4*  Ta- 
khtamirov  and  Volkov41*"52  used  a  simple  model  of  an  IF  to 
relate  the  IF  parameter*  to  IF  grading  and  derived  a  more 
general  form  of  the  IF  potential  within  the  EFA. 

In  a  number  of  papers,  the  theory  of  Refs.  35-39  has  been 
applied  using  the  perturbative  k  ■  F  EFA  approach  to  model 
band  gap  data  and  to  explain  spin  relaxation 
experiments.2"  '31  In  Lau  et  a/.,30  both  the  k  ■  P  term  and  the 
IF  effects  are  treated  as  a  perturbation  on  the  “unperturbed"' 
SL  bands  at  T  and  the  HJI-SO  IF  coupling  (see  Refs.  4.1,  52) 
is  not  included.  The  EFA  k  F  method  in  conjunction  with  the 
IF  Hamiltonian  of  Krebs  and  Vorsin35  (labeled  HBy)  demon¬ 
strated  improved  agreement  wilh  data,12'23 

Wang  et  af.14  and  Magri  ei  rt/.15  19  critically  explored  the 
limits  of  the  EFA  in  specific  applications  by  comparison  to 
the  more  numerically  intensive  empirical  pscudopolentia! 
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(EPS)  Calculations.  For  example,  Magri  et  a7.,g  show  that  the 
standard  EFA  and  EPS  calculations  agree  if  the  EFA  input 
parameters  come  from  the  same  EPS  band  structure  calcula¬ 
tion.  The  electronic  structure  of  InAs/GaSfr  SLTs  lias  also 
been  calculated  via  the  superlattice  empirical  pseudopotem 
rial  method  (SKl'M).^-26  Unlike  the  EPS, 17  the  SEPM  (Ref. 
26)  assumes  that  IP’s  remain  sharp  and  bulk  like,  so  that  the 
method  does  not  differentiate  between  different  types  of 
bonds  possible  at  interfaces.  Nevertheless,  good  agreement  is 
found  with  experimental  absorption  data23"26  on  the  band 
gaps  for  a  series  of  samples,22,23 

In  the  more  accurate  atomic  EFM  (AEPM)  model.,14'19 
the  pseudopotential  is  fitted  to  bulk  bands  and  effective 
masses,  and  the  spin-orbit  interaction  is  included  via  the  non¬ 
local  part  of  the  pseudopotenlial.  An  adjustable  scale  factor 
multiplies  the  kinetic  energy  operator  for  better  fits  to  the 
bulk  band  gap  data  and  to  the  momentum  matrix  coupling 
the  valence  and  conduction  bands.  Without  additional  infor¬ 
mation  on  the  atomic  detail  of  the  IPs,  the  AEPM  (Ref.  14) 
has  limited  success  in  predicting  the  correct  band  gap  trends 
with  layer  thickness.26  Indeed,  AEPM  (Ref.  27)  can  predict 
wrong  dependence  of  InAs/GaSb  band  gaps  on  GaSb  layer 
thickness/*-  However,  Magri  et  at.  found  a  very  good  agree¬ 
ment  with  experimental  band  gaps  for  a  series  of  JnAs/GaSb 
SUs  via  EPM  by  accounting  for  interfacial  diffusion.17 

EPS  and  EFA  represent  two  valid  paths  to  the  goal  of 
modeling  hetcrostnicturcs.  EFA  has  greater  physical  appeal, 
is  vastly  less  computationally  intensive  in  its  setup  and  ex¬ 
ecution,  can  be  more  readily  implemented  by  a  nonspecialist, 
and  is  more  easily  used  for  calculations  involving  electric 
and  magnetic  fields.  Auger  lifetimes,  ionization  coefficients, 
etc.  The  EPS  directly  accounts  for  all  possible  interface 
bonds  and,  given  an  atomic  model  of  interdiffusion,  it  can 
model  interdiffused  IP's.  The  present  EFA  models  these  ef¬ 
fects  through  powerful  symmetry  arguments  about  the  form 
of  the  IF  potential  If  atomic  profiles  are  known  either 
experimentally 54,53  or  theoretically,17  tliey  may  be  included 
in  the  EFA  as  explained  in  Ref.  47  or  Refs,  50-53. 

In  this  paper,  we  investigate  NCA  IF  effects  within  a 
modified  8  X  8  EFA  by  soiving  the  EFA  Hamiltonian  nonper- 
turbativeiy,  based  on  the  solution  of  the  relevant  boundary 
conditions  (BC  approach).  We  adopt  the  IF  Hamiltonian  of 
Ivchenko,  Kaminski,  and  Rdssier4"  and  that  of  Takhtamirov 
and  Volkov,31  including  the  IF  coupling  to  the  SO  band4-10" 
This  paper  is  organized  as  follows.  Section  II  provides  a 
summary  of  existing  approaches  to  the  problem  of  represent¬ 
ing  IPs  in  I  he  EFA.  Next,  in  Sec.  I  ff  the  solution  of  the 
problem  is  formulated  using  the  numerically  stable  represen¬ 
tation  of  the  secular  equation  in  lerms  of  tangents  only, 
which  is  algorithmically  simple  and  guards  against  the  oc¬ 
currence  of  large  exponentials.22  This  method  also  affords  a 
number  of  physical  insights  into  the  influence  of  lFns  on  the 
electronic  structure  of  SL"s  and  the  large  role  of  the  SO  band. 
Wc  find  that  the  IF  coupling  or  the  HH  band  to  the  SO  is 
vEry  important  in  spite  of  the  large  SO  splitting  in  In  As/' 
Ga$b  $U$,  This  is  so  because  the  SO  wave  function  is  con¬ 
centrated  at  the  IP’s  and  is  further  admixed  by  the  LH-SO 

interaction  in  the  k  p  Hamiltonian  at  F,  Further  because  of 
the  singuLar  nature  of  the  perturbation  introduced  by  the  IFns, 


we  show  that  the  problem  cannot  be  dealt  with  perturbativdy 
since  the  solutions  of  the  unperturbed  problem  cannot  repre¬ 
sent  the  rapidly  evanescent  components  of  the  wave  function 
at  the  IFTs.  We  demonstrate  this  by  performing  a  14-band 
perturbative  k  V  EFA  calculation  with  the  IF  interaction  as 
the  perturbation. 

In  Sec.  IV,  we  compare  the  results  of  the  present  calcula¬ 
tion  to  the  band  gap  data  from  several  laboratories,  including 
ours.56  We  explain  the  band  gaps  as  a  function  of  GaSb  and 
InAs  widths  in  terms  of  variations  of  the  heavy -hole  and 
conduction  band  bandwidths.  Wc  show  that  the  cutoff  wave¬ 
lengths  can  be  reduced  by  increasing  the  GaSb  layer  thick¬ 
ness.  Conclusions  are  presented  last  in  Sec.  V. 

IL  REVIEW  OF  INTERFACE  REPRESENTATION  IN  EFA 

For  noncommon-atom  Ill/V  1F*s+  [110]  and  [—110]  direc¬ 
tions  are  inequivalent.  not  being  connected  by  a  fourfold 
rotoin  version  operation  ICAzt  so  that  the  symmetry  at  a 
single  IF  is  reduced  from  D2d  to  C2u  .6,T  As  a  result,  the  z 
component  of  the  total  angular  momentum  J;  is  no  longer  a 
good  quantum  number  so  that  HH  and  LH  arc  coupled  even 
at  f .  To  first  order,  the  5-like  conduction  electrons  aie  unaf¬ 
fected.  The  EFA  representation  of  Dun  deals  with  the  valid¬ 
ity  of  EFA  in  the  IF  regions.'17,4*  I  lore,  three  exiensions  of  die 
theory  to  IF  effects  in  die  EFA  are  described. 

A.  Krebs  and  Vuisin 

Krebs  and  Voisin13  showed  that  the  gianf  in-plane  optical 
anisotropy  of  NCA  group-1  JJ  SL’s  can  be  explained  by  the 
reduction  of  the  assumed  symmetry  lo  the  C2lr  symme¬ 
try.  Using  bond  counting  arguments*  they  showed  that  the 
resulting  symmetry  reduction  couples  the  light  and  heavy 
bands  even  at  f\  unlike  in  the  standard  EFA.32,34  The  effect 
can  be  modeled  in  the  EFA  via  a  short-range  potential  at 
SL  IF  s  42p52  The  short  range  of  the  interaction  in  the  EFA 
was  subsequently  proved  rigorously  by  Foreman47  and  Ta¬ 
khtamirov  and  Volkov.32  However*  the  Krebs- Voisin  theory 
contains  diagonal  HH-llH  and  LH-LH  coupling  with  the 
same  coupling  constant  as  the  HH-LH  coupling.  Such  cou¬ 
pling  does  not  appear  in  the  theory  of  Ivchenko  etal.42'4* 
and,  in  the  theory  of  Takhtamirov  and  Volkov. **' -51  it  ls 
present  only  for  nonabrupt  IF*s  as  an  isotropic  interaction 
with  an  unrelated  coupling  constant.  Recently,  Krebs- Voisin 
theory  was  used  to  calculate  the  hole-spin  relaxation  times 
and  absorption  spectra  in  InAsv'GaSb  SL’s.23~3D  The  strengths 
of  the  anisotropi  c  IF  interactions  (one  for  each  IF  type)  in  the 
theory  were  fit  [o  the  experimental  values  of  the  hole-spin 
relaxation  times  and  band  gap  energies. 

B.  Ivditnko,  Kaminski,  and  Rosskr 

Ivchenko,  Kaminski,  and  Rossler42  represented  (he  effect 
of  IF  anisotropy  on  the  EFA  Hamiltonian  based  on  the  Iheory 
of  invariants  in  which  the  Hamiltonian  is  expanded  in  a  sc¬ 
ries  of  operators  that  are  invariant  under  the  symmetry  op¬ 
erations  of  the  system.^’41  The  value  of  the  expansion  coef¬ 
ficients  is  fixed  by  experiment  but  the  form  of  (be 
Hamiltonian  is  governed  by  the  symmetry  of  the  system.57,53 
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We  adopt  the  form  of  the  IF  Hamiltonian  that  is  given 

by^2 


V=IIxra0S(  r)0, 

(la) 

where  in 

the  |JAfy)  basis 

0 

0 

r/v3 

0 

0 

M\ 

0 

0 

0 

j/V3 

0 

0 

-f/V3 

0 

0 

0 

*> 

0 

0  = 

0 

0 

0 

J2/3 

0 

0 

0 

0 

V*2/3 

0 

0 

0 

0 

0 

0 

a  / 

(lb) 

and  the  j JMj)  basis  is  arranged  in  the  order  |iT&,  tijX  \b 
||  -|)p  H.-i).  au  is  a  lattice  constant  (here,  that 

of  GaSb),  and  HX7  is  the  strength  of  the  IF  potential.  The 
interaction  potential  in  Eq.  (la)  contains  HFbLH  and  HH-SO 
coupling  and  no  diagonal  IIH-HH  or  LH-LH  lerms,  so  that 
the  H&f  description  of  Krebs  and  Vorsin35'30  is  not  equivalent 
to  drat  of  Ivchenko  ct  ai*2  It  is  al*o  common  lo  use  the 
dimensionless  coupling  parameter  r,  defined  through  HXy 
=  hzif2m^al. 

O  Takhtamirov  and  Volkov 

Takhiamirov  and  Volkov51*"5^  adopt  a  particular  model  of 
IF  interaction  in  order  to  derive  an  EFA  Hamiltonian  with  IF 
effects.  In  their  model,  the  EFA  potential  U(r)  across  a 
single  AB  IF  is  graded  as52 

U(t)  =  U4{r)^G(r)[UB(r)-UA(r)l  (2) 

and  all  IP-related  parameters  are  expressed  in  terms  of  the 
Fourier  harmonies  (form  factors)  of  C(r)  and  £/j?(r) 
—  Ua{t)  [see  Eq.  (25)  of  Ref.  52 ],  where  G( r)  is  a  grading 
function  periodic  in  the  in-plane  direction.  As  such*  the 
theory  contains  no  ad  hoc  adjustable  parameters  for  the 
strength  of  the  IF  interaction.  This  theory  naturally  intro¬ 
duces  diagonal  band  coupling  for  nonabmpt  1Fhst  which  be¬ 
come  important  for  valence  and  conduction  bands  in  narrow 
wells.  In  the  limit  of  sharp  IFns+  the  theory  of  Takhtamirov 
and  Volkov  reduces  to  that  of  Ivchenko,  Kaminski,  and 
Rossler,  Eqs.  [la)  and  (lb). 

LIL,  FORMALISM  DEVELOPMENT  OF  THE  MODIFIED 
EFA  THEORY 

This  section  develops  the  mathematical  formalism  for  the 
modified  EFA  theory+  an  EFA  theory  that  includes  IF  effects 
nonperturbatively  through  the  boundary  conditions. 

A.  Hamiltoaiian 

The  notation  used  in  this  paper  is  consistent  with  that  used 
previously/'1  The  super!  atl  ice  consists  of  alternating  layers  of 
material  A  (say,  GaSb)  of  width  2a  and  material  B  (here, 
InAs)  of  width  2b,  In  the  zeroth  unit  cell,  material  A  ranges 


between  and  material  B  between  (™b<2 

**b).  The  period  d^  2a  +  2h. 

In  each  layer  of  a  heterostmeture,  the  ArXjY  EFA  Hamil¬ 
tonian  (here,  N=2X4  =  8)  is  expanded  in  powers  of  kt 
(Refs.  32,  58,  and  59} 

H(kii,ks)=kzH2(k\\)*i+W  3  |(kl)Jtt+/f(,(kl), 

(3) 

where  k7=-  idtdi  and  ffn(k,),  //3(k|}*  and  JT2(Ki)  are 
I  lermitian  and  i  dependent;59  strain  is  included  as  described 
before.52  At  IF *5  there  are  additional  IF-locaiized 

terms  (la), 

V=[HBxAYaiS{z-b)-H^3{2^b)\%.  (4) 

For  InAs/GaSb  £L\s,  there  are  two  possibly  different  lFas — 
the  BA  interface  GaSb-on-lnAs  at  z^b  and  the  interface 
InAs-on-GaSb  at  z=  —  b — either  of  which  can  be  grown 
InSb-tike  Of  GaAs-like.  In  practice,  IP's  often  turn  out  to  be 
mixed  and  interdi (fused,  depending  on  growth  conditions 
and  degree  of  growth  control. S4,5S  For  very  thin  SL's,  a  scalar 
diagonal  IF  interaction  may  become  important*2  and  may  he 
easily  incorporated  in  the  present  theory.60 

Away  from  IP’s  in  each  layer  of  the  £L,  the  EFA  Schro- 
dinger  equation  for  the  envelope  functions  F  (an  N  vector)  is 
given  by52 

d*F  ,  dF 

^  +  i(IF2)-'lIi~iH2r'(H^E)F^  0  (5) 

and  the  solutions  have  the  form 

2N 

k|  *q  -  X  cu  ( ■ kfh  kt)  ik  f*)  ^  j(  k| ,  g  ) ,  (6) 

r  =  Jf 

where  the  N  vector  C[r(  k3  is  the  right  eigenvector1’9  of 

(H2k2  Hq)C=EC  (7) 

at  energy  E  and  exponents  k,  (i  =  ],.„,2JV)  are  the  roots  of 
the  determincntal  equal  ion 

(a) 

In  matrix  notation,  the  column  vector  for  the  envelope  func¬ 
tion,  (6)  is  given  by 

F(z)  =  C  exp(  iKz)c ,  (9) 

where  Kij=kfStj ,  so  that  the  exponential  matrix  is  also  di¬ 
agonal 

[  exp{  iKz)]fj=  exp(  kfz  )StJ.  (10) 

It.  Boundary  conditions 

At  interfaces,  EFA  requires  the  continuity  of  the  envelope 
function  and  current52,55,5*  In  absence  of  IF  effects,  the  two 
continuity  conditions  require  the  continuity  of  the  following 
compact  expression:32 
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f  I  0j|f 
\iNtfl  H2)  \F* 


or,  equivalently,  of 


(I  la) 


M  exp  (iKz)cf  (lib) 

where  the  material-dependent  matrix  M  h  given  by 


for  complex  exponents  (evanescent  solutions),  Instead,  we 
define  A  =  — qdf2)  and  i?  =  tan  and  use  the  gen¬ 

eral  Operator  identity*1 

exp(  2  r'X)  =  ( I  +  J  tan  £)/(  I  -  i  tan  X)  (20) 

and  obtain  two  equivalent  equations  secular  for  matrices  0  , 
and  fl  i 


M— 


(m 


I  .|[C 
i^/2  H2)[CK) 

At  z  =  b,  integrating  the  Hamiltonian  across  the  BA  interface, 
the  two  boundary  conditions  become 


U+i&BA)MA  «P UKAb)cA=MB  zxpUKBb)cB , 


(13a) 


which  can  be  written  as 


[-/(1-(5)e*p(r<&Sj(Kl  -  M)+;(l+ff)«xp{i4ilj) 
X(l+M}]*=  0,^=0,  (21) 

[-i(l  -M)exp(  -  <'$**)(  1  ~iB)+  H I  +  L4)exp(”  !<!»,,,) 

x(i  +  /a)]r=n1y=o>  (22) 

with  the  eigenvectors 

X=  SQ${\A(i-qell2)exp[i{AA - ql)d!2]MAcA ,  (23) 


zx?{i<bSA)MAfi\p{iKjib)cA  =  MgtXf>{iKB!>)c&, 

03b) 

where  the  addition  of  the  ^function  interaction  terms  at  the 
IF  gives  rise  to  a  phase  factor  expfid*^)  with  the  ^phase 
angle"  defined  through  the  interface  matrix 

e  □)■  o+) 

Appendix  A  lists  the  properties  of  the  interface  matrix  and  of 
the  exponential  exp 

At  \htAB  interface  z=  -bT  the  boundaty  conditions  aug¬ 
mented  by  the  Bloch  periodicity  condition,  can  be  written  as 

exp(  -tqd) cxp( i&AB)MA  mp[iXA(d-h)]cA 


y=(co5Aj£)AVj.  (24) 

For  complex  exponents,  Eqs.  (21)  and  (22)  grow  only  as 
hyperbolic  tangents.  For  a  more  symmetric  appearance,  the 
results  are  also  valid  with  the  use  of  A**Xzm{AAa-qct/4)  and 
B  =  tan(Agb—qd/4)„  Note  also  that  the  secular  equation  is 
not  a  simple  function  of  the  sum  or  difference  of  the  FF  terms 
®  ab  and  <PBa  . 

The  energy  eigenvalues  are  at  the  zeroes  of  either  deter¬ 
minant 

|fttl=|ft2l  =  0»  (25) 

since  the  determinants  are  equal.  In  the  absence  of  IF  effects, 
Eq.  (25)  reduces  to  the  simple,  tangents-only  form  found 
earlier'" 


=  MB£xp(-iKBb)c&%  (15) 

where  the  strength  of  the  IF  term  is  §).  For 

symmetric  \F%  JI^=Hyy  in  Eq.  (4). 

a  Secular  equation  and  solution 

The  two  boundary  conditions  can  be  solved  simulta¬ 
neously  to  yield  the  equation 

[exp(i^B)exp[2>(Ajtf-^/2)] 

-  exp[  -  2  i  A  $b  ]  exp(  i  d*  )  ]  V  =  0,  (1 6) 

where  the  eigenvector  can  bo  expressed  as 

F  =  exp(  -  i<bBA)exp(iABb)Mac3  ( I V 


|M  +  5||  =  ||tan(  AAa  —  q  d/4)-\-  tan(  Asb  —  ^cfM)!!  —  0, 

(26) 

After  the  zeroes  of  say,  f]£ i  ,t|.  are  identified,  rhe  correspond¬ 
ing  eigenvector  X can  be  used  directly  to  find  the  wave  func¬ 
tion  expansion  coefficient  cA  via  Eq.  (23);  then,  c3  is  found 
from 


A^[exp(-  z^^)exp(  rAfl&  )  +  exp(  -  i®A3) 

XexjK  -  iApb}]\fBcBf2  (21) 

or,  equivalently,  from 


X=  exp[  -  j  +  2] 


/  -h  sin  |  — j  tan  A $b 


or  as 

r=sxp(lAAb)MAcA;  (IS) 

therefore,  F  can  be  used  to  find  both  cA  and  .  The  A 
matrices  ate  defined  via 


X(MS  cos  KBb)cB.  (28) 

Using  D 2  to  find  T  and  then  cs  h  formally  correct  but 
doubles  the  computational  effort. 

Lastly,  for  symmetric  IP's,  one  form  of 

the  secular  equation  is 


exp(  i Aa  a )  s  ma  cxp(  i  A',,  a }  MA  1 .  (19) 

Although  Eq.  (36)  can  serve  perfectly  well  as  a  secular  equa¬ 
tion,  it  contains  exponentials  that  might  become  very  large 


||exp(j exp(  - i*)H-  B||  =  0.  (29) 

For  completeness,  the  secular  equation  can  also  be  put  in  the 
Kronig-Penney-like  form 
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||[exp( :  $  j  ff)expf  2  i  A  ja)exp(  -  A  Bb) 

+  exp(  -  2  i  A  jA)exp(  i$BA  )cxp( -2l\Aa) 

cos^dJ-O.  (30) 

Equations  in  this  section  can  be  derived  using  the  transfer 
matrix  approach,*1'1  Appendix  B- 

D.  Symmetry 

The  choice  of  origin  in  the  middle  of  layer  B  is  arbitrary, 
Choosing  the  origin  in  the  middle  of  layer  A  (or  translating 
the  system  by  half  the  period),  shows,  that  there  is  a  hidden 
symmetry  i.cM  A  and  B  ate  dummy  labels. 

For  symmetric  IP’s,  Eq.  (29),  changing  the  sign  of  the  inter¬ 
action  at  each  IF,  and  interchanging  the  dummy 

labels  A~B  shows  that  the  determinant  is  independent  of 
the  sign  of  the  interaction.  For  in  equivalent  IF%  Eqs.  (21) 
and  (22)*  if  the  sign  of  the  interaction  at  each  IF  is  changed, 
and  <frw— ►  —  f  and  the  dummy  labels 

A^B  are  interchanged,  it  is  found  that  0,  so  that  the 

eigenvalue  problem  is  again  independent  of  the  sign  of  the 
interaction.  Appendix  C  provides  explicit  solutions  for  a 
number  of  analytic  models  of  SL's  with  symmetric  interface 
interactions;  invariably  the  solutions  are  proportional  to  the 
square  of  the  dimensionless  coupling  parameter  /. 

The  symmetry  group  for  a  symmetric  SL  is  with 
eight  symmetry  operations  of  CV,  -{E.C^JCja 
augmented  with  four  operations  that  interchange  the  two  IP’s 

/Ct,,  (lCtz)~\  Ci])t  Cj Jf.i7  The  factor  group  at  f  is  D14 
itself  with  two  irreducible  representations  compatible  with 
spin  T t  and  FT.  Along  the  q  axis,  the  factor  group  is  C2i} , 
which  has  a  single  two-dimensional  representation  F 5 „  so 
there  is  no  spin  splitting.  Along  the  k T  axis,  the  factor  group 
is  C-2  consisting  of  identity  and  a  twofold  rotation  about  the 
x  axis,  with  two  singly  degenerate  representations  f  j  and 
so  that  all  bands  arc  spin  split.  Representations  F*  and 
F7  at  F  decompose  into  F5  and  F4  along  die  axis,6,7 

For  an  asymmetric  SL,  the  point  group  symmetry  is  Co* 
The  factor  group  at  F  and  along  the  g  axis  is  C2t,  itself  with 
a  single  double  degenerate  representation.  And  along  the  £T 
axis,  the  factor  group  is  C*  with  one  singly  degenerate  rep¬ 
resentation  so  that  the  bands  anticross.  Of  course,  time 
reversal  symmetry  dictates  that  energy  bands  be  degenerate 
upon  reversal  of  wave  vector  and  spin. 


Er  Hermiticlty 

In  Ref.  32,  the  following  identity'  was  proved: 

7AT=A*3  (31) 

Where  the  row  (or  column)  interchange  matrix  is  given  by 


r= 


02) 


which  helped  prove  that  A  +£',  Eq,  (26),  can  be  made  Her- 
mitian  by  interchanging  its  first  and  Last  Ar  rows  (or  col¬ 
umns).  Here,  in  addition,  wc  find  that 


a  L-.D2  D.04  D.DB  0.02  0.1 

4  Y'  \  (Bchr)  + 


FIG.  I .  The  band  structure  of  a  test  IfiAs/GaSb  supcrlartice  with 
44,  \6  A  GaSb/55-46  A  EnAn  without  IF  effects  using  the  standard 
SXS  EfA  (Ref  32).  The  zero  of  energy  is  at  die  CB  edge  of  bulk 
InAs. 


r(/+i*)r=(/-’i*)+  (33a) 

or 

rexp(l*l3)7’=[exp{“  j'd>)]  +  .  (33b) 

The  same  proof  with  the  present  secular  equation,  (21) 
proves  that  H|  cannot  be  made  Hermit ian  because  of  the  IF 
terms,  However, 

m]r=il2+;  (34) 

therefore,  using  Eq.  (25  ^ 

(35) 

so  that  the  determinant  is  real,  an  important  consideration  for 
numerical  purposes.  Lastly,  it  Follows  from  Eq,  (34)  that  if 
O  LX=  AA  and  Q2  then  and 

(TX)+ih  =  \*(TX)+.  (36) 

The  eigenvalue  problem  is  best  solved  by  diagonalizing 
fl  |  as  a  function  of  energy  for  a  given  wave  vector.  The 
determinant  given  by  the  product  of  diagonal  elements*  Aj 
( I  —  l,r„2N),  is  zero  at  an  energy  eigenvalue,  and  the  wave 
vector  Xj  corresponding  to  the  zero  eigenvalue,  A^-G,  is 
used  to  find  the  envelope  function,  Eqs.  (23)  and  (27), 

F.  Typical  hand  structure 

Wc  adopt  a  test  superlattice  for  displaying  typical  results 
of  the  calculation.  The  test  structure  is  an  InAs/GaSb  SL  with 
44.16  A  GaSb/55.46  A  InAs  and  bnth  InSb  IFs.  The  physi¬ 
cal  input  parameters  are  given  in  Table  IX  in  Appendix  D. 
For  future  comparison,  the  band  structure  of  this  SL  without 
IF  effects  (standard  EFA)  is  shown  in  Fig-  1  and  has  the  band 
gap  of  111,4  meV.  The  unstrained  conduction-band-valence- 
band  overlap  is  taken  to  be  0.1 40  eV  for  the  test  case.  Three 
valence  bands  and  one  conduction  band  arc  shown.  The 
bands  labeled  HH  and  LH  do  not  interact  at  T  but  are  hy- 
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bridized  along  the  axis;  HH  bands  are  flat  along  the  q  axis 
because  of  the  small  HH  wave  function  overlap  from  one 
period  to  another. 

Figure  2(a)  shows  the  band  structure  of  the  same  test 
InAs/GaSb  $L  as  in  Fig.  I  but  with  symmetric  IF's  charac¬ 
terized  by  //  vr-  525  meV  or  i  =  5- 1 2>  a  number  on  the  order 
of  magnitude  used  by  previous  The 

resulting  band  gap  is  now  9b. 2 5  meV,  The  bands  along  the  £T 
axis  are  spin  split  and  there  are  a  number  of  avoided  band 
crossings  that  lead  to  level  repulsions,  among  them  a  rise  of 
the  HHl  bandj  thus  a  smaller  band  gap.  Starting  at  T+  the 
bands  are  split  in  opposite  directions,  so  that  the  top  of  the 
valence  band  is  away  from  F ,  Fig.  2(b).  Along  the  q  axis,  the 
bands  are  doubly  degenerate,  and  level  repulsion  results  in 
the  top  of  die  valence  band  being  displaced  from  3'  to  L 
=  {0finTT;d).  Since  the  &-like  Cl ’band  (not  shown)  is  not 
directly  coupled  by  the  perturbation,  its  is  largely  unaffected 
(except  through  indirect  coupling  to  LH  and  SO  bands.) 

Next.  Fig.  2(c)  shows  the  test  SL  with  asymmetric  IP's, 
Mxr=  525  meV  atone  IF  and  ffx r=0  at  the  other.  Along  the 
q  axis  the  hands  are  doubly  degenerate  and  the  Jill  1  band  is 
much  flatter  than  in  Fig.  2(b)  Along  the  axis,  the  bands 
are  spin  split  again  and,  belonging  to  different  symmetry 
representations*  they  all  repel. 

In  Fig.  2(b),  the  positive  curvature  of  the  HHl  band  along 
the  q  axis  requires  explanation.  For  small  q ,  the  nondegen¬ 
erate  perturbation  theory  with  the  perturbation 

\^-^qp]+HXY^[5{z-b)-S{i+b)}&  (37) 


(i.e,*  k  \>  and  IF  terms)  would  require  fourth  order  terms 
q1!2  .  Therefore,  we  use  b)  —  £(z+ £>)]0  it¬ 

self  as  the  perturbation  on  the  already  calculated  SL  bands  as 
a  function  of*/.  Perturbation  V  couples  HH3  (symmetric  with 
respect  to  the  centers  of  GaSb  and  InAs  layers)  and  LH2 
(antisymmetric  about  the  centers  of  InSb  and  GaSb  layers) 
bands  at  the  zone  center.  Near  F,  we  approximate  the  unper¬ 
turbed  HHl  as  (dispcmionlcss)  and  LH2  as  E^_n2 

+  hiq1/2mla2i  where  Fig.  I,  and  find  In  first-order 

non  degenerate  perturbation  theory  that 


£srnt(f  )=  ^hhi  + 


l^ul* 


l^l2  *V 

(3S) 


where  Fsi={HHI  |P|LH2)  and  superscript  zero  denotes 
quantities  evaluated  at  F,  For  symmetric  IFTsh  both  lF?s  con¬ 
tribute  equally  to  the  overlap  integral  From  Eq.  (38)  we 
find  the  following. 

(a)  Since  £hhl>£lhi>  interaction  with  the  LH2  band 
raises  the  FI  Hi  band  in  agreement  with  Fig.  2,  independent 
of  the  sign  of  HXY. 

(b)  HHl  band  acquires  some  of  the  curvature  of  the  LH2 

band  —  336.08  ineV).  Since  the  LH2  band  has  a 

positive  curvature  along  the  q  axis,  the  dispersion  less  HHl 
band  also  acquires  positive  curvature.  Fig.  2(b). 


ct  £Bobf  f1  £B*h*f 1 


ej  fl.02  vM  q-k  c  ns  o.i 
^  £B*h if'  k^  {Bohr)"1 


FIG.  2-  The  energy  band  structure  for  the  lest  Stiperlaltice  using 
the  modifled  5X8  EFA  model  based  on  the  solution  of  the  bound¬ 
ary  value  problem,  Eqs.  (21 )  or  (22).  Along  the  [  1 00]  direction,  rhe 
bauds  are  spin  split.  In  the  [OEM]  direction,  the  bauds  remain  doubly 
degenerate.  From  top  to  bottom,  the  bands  derive  from  HHla  LHI, 
and  HH2  bands  In  Fig.  I  and  aic  strongly  intermixed.  The  irreduc¬ 
ible  representations  are  labeled  for  the  given  directions  using  the 
notation  from  G,  F.  Kosten  J-  O.  DimmocL  R,  G-  Wheeler,  and  H. 
Statz,  Properties  of  the  Thirty  Jlvo  Foml  Groups  (MIT  Press,  Cam- 
bridge*  1963).  (a)  Symmetric  $L  (point  group  f)ld)  characterized  by 
Hxr-  525  meV  or  i-  5. 12  for  both  IP’s;  spin-split  bands  belong  to 
different  representations,  (b)  Valence  band  detail  for  the  symmetric 
SL.  (c)  Asymmetric  SL  (point  group  C^)  cbaiaeterized  by 
=  525  meV  for  one  IF  and  Hxi=0  meV  for  the  other;  spin-split 
bands  belong  to  the  same  representation. 
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TABLE  I.  The  first  few  hands  ^in  meV)  of  die  test  SL  at  F  as  a  funciion  of  die  number  of  bands  included  in  the  k  ■  P  pertuibafive  solution. 
The  second  column  gives  6  of  15  bands  used  as  input  to  the  k  F  perturbation  expansiyn- 


3 

5 

6 

7 

10 

13 

14 

8XS  EFA 

Bands 

Bands 

Bands 

Bands 

Bands 

Bands 

Band:-: 

Baud 

Input 

k  P 

kP 

k  P 

k  F 

k  P 

k  P 

k  P 

LI  13 

-435.46 

-441.06 

-441.06 

-399,57 

-395.01 

-39527 

LH2 

-336.08 

-399.75 

-399.57 

-384.63 

-3S3.80  ‘ 

-382.80 

HH3 

-253.93 

-2S5.I9 

— 193.50 

-193.49 

-193.49 

-190.30 

-191  JO 

HH2 

-40.32 

-67.13 

-64.44 

-64.44 

-55.97 

-55.73 

-55.73 

LHI 

-31.96 

-31.96 

-5.35 

“3.25 

-3.25 

10.43 

10.92 

10.92 

HHI 

94.53 

93.76 

93-77 

97.77 

97.77 

97.77 

98.17 

97.2B 

(c)  At  L ,  HH 1  is  antisymmetric  about  the  center  of  JnSb 
layer  and  LH2  is  symmetric  about  the  center  of  GaSb  layer, 
so  that  the  HH1-LH2  integral  Ft2  is  again  nonzero.  HHI 
hand  is  shifted  upward,  and  the  HHI  hand  acquires  the  nega¬ 
tive  curvature  oF  the  LH2  band  at  L  [compare  Fig.  I  to  Fig. 
2(a)]. 

(d)  On  general  grounds,  anisotropic  IF  interactions  force 
the  band  gap  to  be  indirect,  with  the  valence  band  maximum 
being  either  close  to  F  or  at  Z.  This  effect  can  he  tested  by 
examining  absorption  tails. 

(e)  For  an  asymmetric  SL  of  Fig.  2(c),  one  of  (he  IP's  does 
not  contribute  to  the  interaction  matrix  so  (hat  the  size  of 
effects  seen  in  Fig.  2(a)  is  reduced.  In  Ihe  limit  of  equal  and 
opposite  Uxy  at  the  IF*sn  HHI  interacts  with  LHL  not  LH2- 

(0  Using  the  same  perturbative  calculation,  the  shift  of 
HHI  due  to  its  interaction  with  £02  (second  spin-orbit  band) 
is  proportional  to  f2/A. 

G,  Comparison  to  perturbative  k  P  solutions 

Paralleling  the  calculation  of  Ref.  30,  the  problem  was 
also  solved  pcrturbatively  by  first  order  perturbation  theory, 
with  Eq.  (4)  as  the  perturbation.  Inputs  to  this  calculation  are 
the  8X8  EFA  envelope  functions  al  F, 

is 

F%(  kj  -  0,q  -  0,1)  =  2 ,  Cp(  k,  =  0*i)  exp(  i kp  > 

Xcf'(k1l=0<?  =  0)  (39) 


and  energies  £^(^=0^  =  0)  for  band  N  from  the  standard 
8x8  EFA.  Using  the  wave  function  expansion  in  layer  B 
(here.  In  As),  the  first  order  perturbation  matrix  between 
bands  A/  and  M  is  given  by 

ritM-aol  vf + O 

X &afl(  1  +  i  tan  kjb  tan  *?/>)+■ i(  fifr-  NBXAY) 
X&Bft(iaiikjb—  tan  k*b)},  (40) 

where 

fi^cotKBb)cB.  (41) 

To  first  order,  the  perturbation  matrix  depends  on  the  sum 
and  difference  of  the  perturbing  IF  terms. 

The  perturbative  solution  was  tested  as  a  fund  ion  of  the 
number  of  bands  included  in  the  perturbative  treatment, 
Table  L  Using  the  test  SL  with  symmetric  IF’s,  up  to  15 
bauds  were  calculated  at  F,  including  Cl  through  C2h  HHI 
through  RH7,  and  LH1  through  LH6,  spanning  the  interval 
from  about  -700  to  600  meV.  The  size  of  the  perturbation 
matrix  elements  |  converged  rapidly,  im¬ 

plying  (hat  (he  degenerate  perturbation  theory  should  also  be 
convergent. 

From  Table  I,  about  14  (28  with  spin)  bands  are  required 
to  converge  the  first  few  hands.  The  HHI  band  shifts  by  only 
3  meV  as  (he  result  of  interaction  with  14  bands,  indicating 
weak  coupling.  The  LH3  band  ls  greatly  affected,  moving  by 


TABLE  II-  First  few  energy  Levels  for  the  test  SL  with  symmetric  IFs  characterized  by  Hxr=525  mtV  calculated  via  the  14  x  (4 
perturbative  k  P  scheme,  Eq.  (40),  and  with  the  present  BC  approach,  Eq.  [29).  Separate  columns  are  for  results  without  IF  coupling  to  the 
SO  or  LH  bands. 


Level 

//vr=o 
&X8  EFA 

k-F 

No  SO 

kP 

No  LH 

k'F 

All 

BC 

No  SO 

BC 

No  LH 

BC 

Alt 

Cl 

205.92 

212.54 

2071 8 

208.18 

216.55 

208.36 

212.23 

HHI 

94.53 

93-42 

94.64 

97.28 

93.47 

100,54 

315.96 

LHI 

-31.96 

-4.E4 

-27.01 

10.92 

-3.91 

-13-09 

57.32 

11 H  2 

-40.32 

-58.80 

-35,87 

“55.73 

-58.68 

-33.16 

-47,64 

11113 

-253.93 

-22*.  19 

-23S.33 

-191.30 

-224.24 

-192.40 

-102,28 

LH2 

-336D8 

-362.84 

—  352.18 

-383,80 

-355.51 

— 345.02 

-367.S5 
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proach  is  less  converged,  There  are  much  greater  differences 
between  the  k>P  and  BC  mclhods  in  the  “no  £”  columns, 
where  only  the  HJJ-SO  IF  coupling  is  retained.  With  all  in¬ 
teractions  included,  the  differences  between  the  k  P  and  BC 
columns  are  very  large.  We  contend  that  the  perturbative  klJ 
method  cannot  adequately  represent  the  LF -peaked  SO  wave 
function,  whose  large  curvature  adds  more  kinetic  energy  to 
energy  eigenvalues  than  is  reflected  in  the  k-F  approach. 

Figure  3  plots  C I ,  HI-31,  Lilt,  and  HH2  eigenvalues  in 
the  14X  14  k  ■  Is  and  8X8  BC  approaches  as  a  function  of  the 
strength  of  the  IF  interaction,  HHS  rises  in  energy  with  in¬ 
creasing  interaction  strength  as  the  result  of  its  hybridization 
with  lower  bands.  The  greatest  difference  between  the  two 
calculations  is  for  the  LH1  level  because  it  contains  the  larg¬ 
est  SO  admixture  due  to  the  direct  LH-SO  coupling  in  the 
k-p  Hamiltonian.  The  HHI  difference  is  22  meV,  which  al¬ 
though  small  on  the  scale  of  the  graph  (or  on  the  scale  of  the 
VB  offset,  i.e.,  depth  of  the  VB  QW),  is  very  large  on  the 
scale  of  the  band  gap  of  1 00  meV,  The  large  difference  in  the 
LHI  energies  would  result  in  greatly  different  Auger  life¬ 
times  predicted  by  the  two  calculations  based  on  exactly  the 
same  Hamiltonian.  Lastly,  calculated  interband  absorption 
peaks-  in  the  BC  approach  would  be  closer  spaced  than  in  the 
It!1  approach. 

Another  way  of  studying  the  influence  of  the  spin-orbit 
band  is  to  consider  Table  111,  which  breaks  down  various 
components  of  the  calculation,  Without  IF  coupling,  H^y 
—  0,  in  going  from  the  4x4  to  6  x  6  and  then  8  x  g.  EFA 
results,  the  very  large  effect  of  the  k-p  coupling  of  the  LHI 
band  to  the  SO  and  C  bands  is  apparent.  With  the  hill  inter¬ 
action  turned  on,  the  34X14  model  gives  lower  energies,  as 
noted  in  Fig.  3.  By  comparing  ibe  three  k  F  calculations,  it  is 
clear  that  the  full  effect  of  the  IF  interaction  is  turned  on  only 
when  the  SO  band  is  included-  With  IIH-IJI  IF  coupling 
only,  the  positions  of  the  three  (and  most  other}  bands  in  the 
BC  8X8  EFA  and  the  14x  14  k-P  perturbation  model  are 
very  similar,  indicating  the  crucial  importance  of  the  IF  cou¬ 
pling  to  the  SO  band. 


TABLE  IB.  The  calculated  energies  (in  meV)  for  (be  test  SL  fit  T  using  several  models  for  the  modified 
EFA  (BC  columns)  and  the  NX  14  perturbative  k-F. 


Interaction 

Level 

4X4  HI 

BC 

6X6  HLS 

BC 

8*8  CHLS 

BC 

14x  E4  hi 
Pertu  ft  [stive 

HJr  13 

94.53 

94.53 

94.53 

94.53 

LHI 

-29.54 

-10.75 

-31.96 

-31.96 

HH2 

-■40.32 

-40,32 

-40.32 

-40.32 

HXif=252  meV 

HHI 

96.16 

320-66 

115.96 

97.28 

All 

LHI 

8.85 

69.62 

57.32 

10.92 

HH2 

-64.46 

-36.76 

-47.65 

-55.73 

Hxr=525  meV 

HHI 

96.76 

97.98 

93.48 

93.42 

Itll-LH  IF 

LHI 

8.85 

22.74 

-3.91 

-4.83 

Coupling 

HH2 

-64.46 

-55.45 

-58.68 

-58.80 

Hxr^525  meV 

HHI 

102.49 

100.54 

9L64 

HH-SG  IF 

LHI 

-8-71 

-1309 

-27.01 

Coupling 

HH2 

-11.66 

-33.16 

—  35.87 

c  aw  zw-  sw  frw 


ii^IikV) 

FIG.  3.  Plots  of  the  Cl.HHL  LHI,  and  III  12  eigenvalues  for  the 
perturbative-  I4x  J4  k  F,  Eq.  (41),  and  the  modified  8  X  8  EFA,  Fq. 
(21h  approaches  as  a  Function  of  the  strength  of  the  IF  interaction 
for  tEie  test  super  lattice 

some  40  meV,  since  symmetric  IP's  couple  opposite  purity 
LI  1-1 1 H  pairs  such  as  LHI  and  FIH2.  Similarly,  HHI  and 
LH2  bands  arc  expected  to  interact,  but  tor  the  parameters  of 
die  test  SL,  these  bands  are  already  separated  by  450  meV. 
The  symmetric  IFhs  have  almost  no  effect  on  the  conduction 
band  (not  shown  in  Table  1), 

Table  II  provides  a  comparison  of  the  calculated  energy 
levels  for  the  tesl  SL,  Figs.  1  and  2,  using  the  14X  14  per¬ 
turbative  k-P  solution  and  the  BC  method  with  the  modified 
8X8  EFA.  For  the  perturbative  solution,  we  use  14  doubly 
degenerate  bands  calculated  at  T — one  conduction  band  and 
thirteen  valence  bands.  One  finds  that  as  long  as  the  SO  band 
is  not  included  (bolh  lino  SO”  columns),  the  results  of  the 
perturbative  k-F  and  the  BC  methods  are  comparable,  al¬ 
though  the  BC  approach  gives  slightly  more  positive  ener¬ 
gies  especially  for  the  deeper  bands  for  which  the  k-P  ap¬ 
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Lastly,  Table-  III  shows  that  the  effects  of  the  HH-LH  and 
HH-SO  IF  couplings  are  not  simply  Additive— the  problem  is 
highly  nonlinear.  An  additional  reason  for  studying  the  4 
X4  and  6X6  models  here  is  that  they  are  devoid  of  possible 
spurious  solutions  that  sometimes  appear  in  the  S  x  8  EFA 
model. 

The  BC  approach  .solves  the  problem  exactly  and  should 
be  used  as  the  standard  of  comparison.  The  BC  method  is 
“horizontal" — each  wave  function  is  expanded  in  terms  of 
properly  joined  solutions  of  the  bulk  problem  at  the  energy 
and  wave  vector  of  interest,  Eq,  (6).  The  perturbation  ap¬ 
proach  is  “vertical" — each  wave  function  is  expanded  in 
terms  of  the  eigenstates  of  the  unperturbed  Hamiltonian  (39) 
that  are  ascending  in  energy.  For  example,,  in  the  BC  method, 
the  wave  function  of  the  HHI  state  has  5>G  components  at 
the  energy  of  the  HHI  state,  components  that  are  very 
sharply  peaked  at  the  IP's,  since  the  SO  state  is  highly  eva¬ 
nescent  at  die  energy  of  die  HHI  state.  Already  in  the  tight- 
binding  calculation  of  Schulman  and  Chang7  (see  their  Table 
VI).  HI!  states  are  shown  to  have  large  LH  and  SO  admix¬ 
tures  at  r. 

Additionally,  in  Appendix  E,  a  test  case  demonstrates  that 
HHI  in  A  model  SL  shifts  more  when  computed  in  the  6 
X6  HLS  model  than  in  the  4X4  HS  and  4x4  ML  models 
combined.  This  implies  that  the  problem  is  highly  nonlinear, 
thus  inappropriate  for  perturbative  treAl merit.  Then,  in  Ap¬ 
pendix  F,  we  prove  that  the  anisotropic  interaction  at  a  single 
IF  is  so  singular  that  it  can  bind  A  hole. 

Overall,  the  present  result  indicate  that  the  spin-orb  il 
band  h  very  important  bolh  qualitatively  and  quantitatively. 
In  particular,  the  position  of  the  HHI  valence  band  is  seen  to 
vary  greatly  when  the  IF  coupling  to  the  SO  band  is  in¬ 
cluded.  Already,  Ivchenko,  Kaminski,  and  Rossler42  showed 
that  the  SO  band  is  important  in  fitting  the  anisotropic  ex¬ 
change  splitting  of  exci  tonic-  levels  in  GaAs/AlAs  SL's  by 
demonstrating  that  the  fit  changes  from  /  =1.4  to  t=  0.6 
without  and  with  the  SO  band  included,  respectively42 

IV.  COMPARISON  TO  EXPERIMENT 

In  order  to  compare  experiment  to  theory,  the  modified 
EFA  requires  parameters  for  the  GaAs-like  and  InSb- 
like  IP’s.  For  a  symmetric  structure,  however,  only  one  Hxy 
parameter  is  required.  For  the  determination  of  this  param¬ 
eter  we  grew  and  measured  the  photoresponse  of  two  sets  of 
InAs/GaSb  SL’s  with  InSb  IFLs — one  with  variable  GaSb 
layer  widths  and  the  other  with  variable  InAs  layer  widths,3* 

A.  Ei  pc  rim  filial 

The  InAs/GaSb  SL's  were  grown  by  molecular  beam  ep¬ 
itaxy  (MBE)  using  elemental  metals  for  Ga,  In,  and  valved 
cracker  cells  for  Sb  and  As.  The  substrates  were  epiready 
{ 1 OO)  Te-doped  GaSb  wafers.  First*  GaSb  buffer  Sayers  were 
grown  on  the  substrates  at  the  growth  temperature  of  480  DC. 
Next,  the  temperature  was  lowered  to  4  1G±  5  °C  to  grow  the 
SL  structure.  Here,  growth  temperatures  are  pyrometer  read¬ 
ings  referenced  to  the  GaSb  oxide  desorption  temperature  of 
530  aC.  The  InSb-like  IF  bonds  were  inserted  between  the 
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TABLE  IV.  The  positions  of  the  first  heavy  hole  (HHI)  and 
conduction  (Cl)  bands  for  two  values  of  unstrained  InAs 
conduction/GaSb  valence  band  overlaps 


JnAs/GaSb 
Overlap  (eV) 

InAs 

(A) 

fiaSb 

(A) 

HHI 

(meV) 

Cl 

(mcV) 

0.140 

S3 

£5 

124.95 

129.99 

86 

125.26 

127.94 

87 

125.85 

123.96 

0.150 

S3 

B3 

334.17 

135.8,3 

84 

£4 

334.69 

153.73 

85 

£5 

134.92 

131.61 

layers  to  reduce  the  SL  strain.  The  number  of  SL  periods  was 
fixed  at  40  and  the  growth  rates  of  GaSb,  In  As,  and  In-Sb 
were  1 .32,  0.38,  0.47  A/s,  respectively.  The  samples  in  the 
same  series  were  grown  consecutively  within  a  short  time 
period  to  minimize  the  variation  of  growth  rales.  The  shutter 
sequence  process  wfithin  the  series  was  kept  the  same  to  en¬ 
sure  consistent  intermixing  between  the  layers.  The  V/1TI 
beam  equivalent  pressure  flux  ratio  was  approximately  3.5 
for  the  In  As  and  2-5  for  the  GaSb  growth.  Cracker  tempera¬ 
tures  w‘ere  900  X’  for  As  and  £50  DC  for  Sb.  These  growl  h 
conditions  were  optimized  for  the  smoothest  surface  mor¬ 
phology  and  the  best  structural  quality.  With  optimized  con- 
ditinns,  the  MBE  growth  process  maintained  the  intended 
period  within  1  A. 

To  determine  the  actual  SL  period  and  strain,  high- 
resolution  x-ray  rocking  curve  (HRXRD)  measurements 
were  performed.  The  actual  periods  were  determined  from 
the  SL  satellite  peak  splitting  in  the  (004)  spectra.  For  the 
present  study,  individual  layer  thicknesses  were  only  esti¬ 
mated  but  arc  expected  to  be  dose  to  the  nominal  values.  It 
is  reasonable  to  assume  that  the  precisely  calibrated  growth 
rates  and  a  consistent  shutter  sequence  process  within  (he 
scries  should  provide  highly  accurate  individual  layer  thick¬ 
nesses  with  consistent  intermixing  between  the  layers. 

To  deientiinc  the  band  gap  energies  and  other  band  struc¬ 
ture  parameters,  Fourier  transform  infrared  (FT1R)  photore¬ 
sponse  spectroscopy  was  used-  The  spectral  response  for  this 
series  was  measured  using  a  BIO-RAD  6000  FTIR  spec¬ 
trometer  system.  The  samples  were  contacted  with  indium 
stripes  and  then  illuminated  from  the  front  at  normal  inci¬ 
dence.  The  speclrA  were  collected  at  10  K. 

B.  Fitting  VB-CB  overlap 

The  EEA  requires  A  few  parameters  to  be  fixed  by  experi¬ 
mental  data.  For  individual  layers,  the  required  data  comes 
from  measurements  on  the  constituent  bulk  components  of 
the  SL,  Table  IX,  Appendix  D.  This  leaves  only  information 
on  the  IF among  them  the  conduction- valence  band  over¬ 
lap.  While  this  offset  has  been  measured  by  independent  ex¬ 
periments  and  calculated  theoretically,  here  the  offset  is  fixed 
by  using  the  fact  that  for  equal  widlh  SL's,  the  metal- 
insulator  transition  takes  place  at  about  S3.5A,^  at  which 
anisotropic  IF  effects  are  negligible.  Table  IV  show's  that  the 
V6-CB  overlap  of  0.150  eV  gives  the  best  fit  to  the  S- M 
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TABLE  V.  Band  extrema  and  the  conduction  band  bandwidth  for  a  series  of  grown  5LFs  with  the  In  As 
width  of  20.5  A  and  the  GaSb  width  ranging  from  IS  to  27  A  and  purposefully  gmwn  InSb  IP's,  The 
unstrained  VE-CB  overlap  of  150  mcV  was  used. 


GaSb  width 

(A) 

HHI 

(meV) 

Cl 

r-r 

4 

(meV) 

Exp. 

(mcV) 

Bandwith 

L-f 

Centroid 

a+nn 

27 

I3L6 

438,1— >626 J 

306.5 

299 

E87.9 

532.1 

24 

129.3 

421.8^641.9 

292.5 

293 

220.1 

531.2 

21 

126.4 

402,3— >660.  B 

276.1 

279 

258,5 

53L  5 

18 

121  8 

378.7— >683.3 

256.9 

254.5 

304.6 

531,0 

transition  point  whereas  the  offsets  of  0. 140  produce  the  S-M 
transition  around  86.5  A.  The  10  sneV  difference  for  the  1 
ML  difference  in  the  layer  widths  gives  an  indication  of  the 
overall  accuracy  of  the  calculation. 

C.  Variable  GaSb  study 

Four  samples  with  GaSb  widths  ranging  from  1 8  to  27  A 
and  the  In  As  width  held  constant  at  2Q.5  A  were  grown  with 
InSb  IFhs  as  described  above.  The  HXy  parameter  character¬ 
izing  FnSb  IPs  was  calculated  by  fitting  the  band  gap  of  the 
18  A  sample  and  then  used  in  the  calculations  for  the  rest  of 
the  SL's  in  this  paper.  The  value  of  Hxy=  580  meV*  or  t 
=  5.66,  fit  the  gap  to  within  a  couple  of  meV.  For  compari¬ 
son,  Lau  and  Flatte30  use  symmetric  and  antisymmetric  in¬ 
teraction  potentials  V$=  500  meV  and  F*=36G  meV,  re¬ 
spectively,  and  Olesberg  er  fit  the  spin  relaxation  times 
in  short-period  InAs/GaSb  SLas  with  the  VA  between  300  and 
500.  In  another  system*  Ivchenko  et  alAl  use  r  =  0.9  and  0.32 
for  AIAs/GaAs  and  obtain  a  tight-binding  estimate  of  / 
=  0.76.  Also*  for  an  InAs/GalnP  IF,  Ivchenko,  Toporov,  and 
Vbisin44  find  —  4.47.  Therefore,  the  present  estimate  is 
consistent  with  the  order  of  magnitude  found  by  others.  The 
results  lor  the  other  SLss  in  the  series  and  the  experimental 
hand  gaps  arc  given  in  Table  V. 

The  agreement  between  experiment  and  theory  in  Table  V 
for  (he  band  gaps  Is  very  good,  with  small  differences  being 
within  the  uncertainties  in  the  layer  widlhs.  The  present 
modified  Sxg  EFA  calculation  shows  that  heavy  holes  are 
largely  confined  in  the  GaSb  layers  while  electron  wave 
functions  overlap  considerably  from  one  InAs  layer  to  an¬ 
other.  As  GaSb  layers  become  narrower*  the  HH-derivcd  top 
of  the  valence  band  decreases  in  energy.;  at  the  same  time, 
the  overlap  between  the  electron  wave  functions  Increases, 


leading  to  an  increase  in  the  GB  bandwidth.  Table  V.  With 
increasing  CB  bandwidth*  the  boltom  of  the  CB  decreases  in 
energy  at  a  rate  greater  than  the  decrease  of  the  top  of  the 
VB*  leading  to  the  observed  decrease  in  the  band  gap  with 
decreasing  GaSb  layer  width.  Also*  the  center  of  the  CB 
remains  nearly  constant,  and  only  its  top  and  bottom  move* 
increasing  the  bandwidth.  This  trend  is  naturally  explained  in 
the  tight-binding  picture  of  band  formation.  The  net  effect;  is 
that  the  bottom  of  the  CB  moves  down  faster  than  does  the 
HH  band*  thus  the  band  gap  narrows. 

The  fact  that  the  CB  rises/falls  faster  than  the  valence 
band  with  changes  in  GaSb  thickness  agrees  with  the  find¬ 
ings  of  Ongstad  ei  al 22  We  supplement  these  findings  with 
our  observation  that  the  middle  of  the  CB  is  largely 
invariant — only  its  bandwidth  changes. 

IX  Variable  In  As  width  study 

Table  VI  presents  the  results  of  the  variable  In  As  width 
study,  with  the  inAs  layer  width  ranging  from  16  to  23.5  A 
and  the  GaSb  width  held  constant  at  24  A.  The  calculation 
uses  the  same  asymmetry  parameter  Hxy=5#0  meV  for  both 
IF’s  as  that  used  in  Table  V  for  the  variable  GaSb  width 
study. 

In  Table  V[,  the  HHI  level  is  largely  invariant  to  InAs 
width  variations  especially  lor  thicker  InAs  layer  widths.  Cl 
varies  in  expected  fashion*  decreasing  rapidly  for  larger 
InAs  layer  widths,  because  of  the  small  InAs  CB  mass. 
The  experimental  result  for  the  16  A  sample  appears  incon¬ 
sistent  with  the  results  for  the  other  four  samples.  The  overall 
agreement  Is  good  with  only  one  adjustable  parameter  ob¬ 
tained  on  a  sample  from  a  different  set.  Table  V.  The  HHI- 
LU I  separation  is  much  smaller  than  the  gap  so  that  Auger 
recombination  cannot  be  suppressed  through  final  state  opti¬ 
misation. 


TABUS  VI,  Calculated  and  measured  band  gaps  for  several  InAs-GaSb  SL's  with  a  constant  GaSb  width 
of  24  A  ajid  the  InAs  layer  widths  ranging  from  16  to  23.5  A.  The  same  asymmetry  parameter  Hxy 
=■580  mcV  as  in  Tabic  V  was  used  for  both  IF's.  The  unstrained  VB-CB  overlap  of  150  mcV  was  used. 


inAs  width 

(A) 

LH1 

(meV) 

HHI 

(meV) 

Cl 

(meV) 

CmeV) 

Exp. 

(meV) 

16 

74.34 

132.67 

503.20 

370.53 

343 

17.5 

70.37 

130.99 

477.67 

346.68 

340 

19 

66.69 

129.71 

454,05 

324.34 

314 

20.5 

63.26 

128.73 

432.16 

303.43 

298 

23.5 

5B.05 

127.36 

392.93 

26557 

269.6 
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E«  19  pm  H4L xct|7 It,  Wei  al  (Ref.  11) 

This  sample — I  ML  InSb  IP's,  17  ML  lnAsh  and  7  ML 
GaSb  -was  grown  with  both  InSb  IP's  anti  is  thus  suitable 
Ibr  simulation  using  the  previously  found  //jr— SEGmeV, 
The  results  here  are  Cl  =  184.13  meV  and  HH3  of  SI 9.63 
meVj  for  the  gap  of  64.5  mcV,  or  a  cutoff  of  19.2  pm. 
The  experimental  numbers  arc  50%  cutoff  at  18.5  ^tm 
(66  mcV)  and  a  0%  cutoff  of  about  20.5  pm  (or  60.5  mcV). 
Also,  LHI  =25.48  meV,  so  HHl-LHl  separation  is  larger 
than  the  gap.  The  standard  8X8  EFA  predicts  a  7,5 
cutoff. 

F.  32  /im  sample,  Wei  ei  ui  (Ref,  11) 

For  this  sample— i  ML  InSb  IFs  16  ML  lnAs+  and 
4  ML  GaSb — we  again  use  i/^=580meV  and  find  Cl  at 
154.38  and  Mill  at  102.43  for  the  gap  of  51.95  meV.  The 
experimental  fit  to  tf0A  (resistance-area  product)  data  vs 
temperature  tor  this  sample  yields  the  activation  energy 
of  45.5  meV.  The  standard  8X8  EFA  cude12  predicts  a 
5,5  jun  cutoff.  Overall,  the  results  for  SL’s  with  controlled 
InSb  IPs  agree  well  with  experiment  both  in  the 
MWIR  and  LWIR  windows  for  samples  where  IP  effects 
should  be  large  owing  to  the  thinness  of  at  least  one  of  the 
layers. 

Gr  Mixed  inkrbai; 

Ongstad  etal.22  and  Kaspi  et  vt23  grew  two  sets  of 
samples  with  mixed  JF*S-  Using  Hx y  =  5 80  meV  for  InSb- 
like  IF*s>  we  use  the  band  gap  of  the  6  ML/6  ML  sample  to 
fit  the  coupling  parameter  for  GaAs-like  IP's,  It"  one  were  to 
fit  both  parameters  at  the  same  lime,  the  graph  of  MlXyb  vs 
would  yield  an  oval  curve*  thus  the  need  to  fix  one  of 
the  parameters  beforehand.  Even  with  one  HXy  fixed,  the 
calculated  band  gaps  arc  A  rapidly  varying  function  of  the 
other  Hxv- 

Table  VII  presents  the  experimental-theoretical  compari¬ 
son  for  the  mixed-IF  samples  where  the  fit  As  width  is  held 
constant  at  6  ML,  The  measured  band  gap  for  the  6  ML/6 
ML  SL  was  well  fit  with  meV  and  H}Xyb  held  fixed 

at  5  SO  mcV,  This  does  not  mean  that  GaAs  JF^s  do  not  in¬ 
troduce  a  perturbation,  only  that  reasonably  good  fits  can  be 
obtained  for  HXy  around  zero  to  within  (he  accuracy  of  the 
calculation  since  A  small  variation  of//^b  can  lead  to  a  large 
variation  of  (c.g.,  Fig.  3).  In  addition,  we  hAve  no 

direct  compArison  between  the  IF**  grown  by  Ongstad  ei  af. 
and  Kaspi  et  a/.  and  those  grown  in  our  lab.  The  small  value 
of  indicates  that  IF  effects  for  mixed  interface  Si's  are 
smaller  that  for  SL's  with  both  InSb-like  TF+s. 

i .  Table  VISE  presents  the  experimeiitaLrheoretieal 
comp arlson  for  the  mixed- IF  samples  of  Kaspi  ei  a!-,  where 
the  In  As  width  is  held  constant  at  8  ML  and  the  theoretical 
IF  parameters  are  held  at  H]£yb=  580  meV  and 
-OmeV  as  in  Table  VIT.  The  true  valence  band  maximum 
for  the  cases  of  Tables  VII  and  VIII  is  slightly  away  from  the 
zone  center  and  not  more  than  0.5  meV  greater  than  the 
center  zone  values  listed  in  the  tables.  In  both  tables h  the 
agreement  with  experiment  is  very  good  and  comparable  to 


PHYSICAL  REVIEW  E  69*  E  5532 1  (2004) 

TABLE  VI I.  Experimental-theoretical  comparison  for  the 
nruxccM  F  series  with  the  a  6  ML  In  As.  1  lere  HXy  =  5  SO  meV  for  Lhe 
InSb  IF  arid  for  the  GaAs  IF,  Unstrained  In. As  CB/ 

GaSb  VR  overlap  of  0.150  eV.  Tlie  band  gap  data  arc  from  Refs.  22 
and  23. 


M  L’s 

IjiAs 

GaSb 

HHJ 

(mcV) 

Cl 

(mcV) 

Et 

ImeV) 

Exp 

(mcV) 

103.7K 

43  8.75 

315-1 

m  PL 

6/6 

18.18 

18.29 

326  Abs 

121.74 

487.30 

365.4 

355  PL 

6/9 

18.18 

27,43 

358  Abs 

130.59 

527.06 

396.5 

3E2  PL 

6732 

18.  IB 

36.57 

386  Abs 

338.81 

566-66 

427,9 

415  PL 

6/IB 

18.18 

54.56 

416  Abs 

142.65 

533.31 

440.7 

425  PL 

6/24 

18,18 

73.14 

427  Abs 

that  found  in  the  SEPM  (Ref  22)  And  FFM  calculations. 1  ’ 
The  results  of  variable  GaSb  width  studies  in  Tables  V?  VII, 
and  VIII  are  presented  in  Fig.  4- 

V,  CONCLUSIONS 

The  band  structure  of  InAs/GaSb  £L*s  was  investigated 
within  the  framework  of  the  8X8  EFA  method,  modified  by 
the  inclusion  of  IF  effects,  based  on  the  solution  of  the  asso¬ 
ciated  boundary- value  problem  And  the  Ivchenko-Kanninski- 
Rdssler  representation  of  the  IF  Hamiltonian.  It  was  found 
that  a  consistent  application  of  the  theory  requires  the  inclu¬ 
sion  of  the  SO  band  and  that  its  inclusion  dramatically  alters 
the  bands  from  those  calculated  with  the  HH-LH  IF  coupling 
alone.  With  the  inclusion  of  the  //-SO  IF  coupling,, 
perturbative  solutions  converge  to  incorrect  limits  because 
(he  SO  admixture  is  highly  peaked  at  the  lF*s  and  the 
IF  interaction  is  singular  enough  to  potentially  bind  a  hole 
at  the  IF. 

For  a  test  InAs/GaSb  SL  with  symmetric  IF’s  and  for 
reasonable  strengths  of  Ehe  IF  interaction,  the  band  structure 
exhibited  large  splittings  and  avoided  crossings.  Analytic  so¬ 
lution  of  the  EFA  model  in  several  important  limits  provided 
additional  physical  insight  into  the  band  formation  in  $L+s 
and  QW's  in  the  presence  of  IF  effects.  The  calculated  results 
were  also  verified  against  these  analytical  models. 

Very  good  agreement  was  found  between  experimental 
results  and  theory  on  several  sets  of  SL's  (both  MWIR  and 
LWIR}  with  symmetric  InSb  IP's,  using  one  adjustable  pa- 
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TABLE  VIII.  Sam*  as  Table  VII  for  the  mixed,- IF  series  with  die 
In  As  laycnsidth  of  8  ML. 


ML’s 

JnAs 

GaSb 

Hin 

(mcV) 

Cl 

(mcV) 

Eg 

(meV) 

Exp 

(mcV) 

111. 74 

385.07 

273.33 

277  PL 

B/B 

24.23 

24.38 

283  Abs 

127.04 

437.60 

33  0.65 

304  PL 

\m 

24.23 

36-57 

308  Abs 

134,38 

464.70 

330.32 

330  PL 

\m 

24.23 

48.76 

332  Abs 

141.17 

486.92 

345.75 

346  PL 

24.-8 

2423 

73  J  4 

351  Abs 

144.20 

493.56 

339.36 

346  PL 

32/8 

24_23 

97.53 

350  Abs 

145.86 

495.57 

349.7 

352  FL 

4m 

24_23 

121.91 

353  Ah* 

rameter  for  the  strength  of  the  IF  interaction  at  InSMike 
IF’s,  Then,  with  one  more  adjustable  parameter  for  CxaAs 
IF'S,  very  good  agreement  was  found  for  two  series  of 
samples  with  mixed  OaAs-like  and  InSb-llke  IF's. 

Overall,  it  was  shown  that  a  consistent  application  of  the 
EFA  method  with  the  inclusion  of  well  established  IF  effects 
can  provide  useful  physical  insights  and  possesses  good  pre¬ 
dictive  capacity  in  the  design  of  NCA  $L’s-  The  formalism 
derived  here  is  algorithmically  simple,  numerically  stable, 
computationally  fast,  and  capable  of  providing  rich  physical 
insight  into  band  formation  in  $L’sa 
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APPENDIX  A;  INTERFACE  MATRIX 

The  2NX2N  interface  matrix  [Fq.  (14)]  was  defined  as 
/  0  G'\ 

<*>  =  ffXy«a\e  Qj.  (Al) 


FIG.  4.  Collected  results  of  the  variable  GaSb  width  studies 
listed  in  Tables  V,  VII,  and  Vlll,  where  the  solid  curves  are  calcu¬ 
lated  with  rhe  modified  EFA  and  the  experimental  data  are  from 
Haugan  er  at  (Ref.  56),  Ongslad  et  at  (Ref  21),  and  Kaspi  et  at 
(Ref.  23). 


where  the  NXN  matrix  0  is  Mermitian.  Matrix  cannot  be 
diagonalized  by  a  similarity  transformation  but,  nevertheless, 
has  several  interesting  properties. 

(1)  Multiplicative  inverse 


(/-»)-= 


/  0 

l-Wj / 


=(/+id>) 


-I 


(A2) 


(2)  Multiplicative  associativity 

(/+^l)(/+i*l)-[/+f(^+*a)J  (A3) 


All  these  properties  are  shared  with  the  exponential  ma¬ 


trix  defined  through  the  Taylor  expansion 

cxp(^)=/-Fi^^(/4>)2/2  +  — (A4) 

since  all  powers  of<b,  Eq.  (Al)p  higher  than  the  first  are  zero. 
Analogously  So  Eqs.  (A2)  and  (A3)>  il  follows  that 

exp(  - Bt>)  =  [eKp{  „  (A5) 

exp(i<f*  |  JfficpO'd*;)—  exp[i(d>]  +d>|)].  (A  6) 

Also*  the  determinant 

|«p(tt*)|=l.  (A  7) 


Moreover*  trigonometric  functions  defined  through  the  expo¬ 
nential  satisfy 


/  o 

o\ 

sind»  =  iand>=  |  ^ 

or 

(A8) 

ccb4>=/3 

(A0) 

and  such  identities  as 


sin2  *b  +  cos2  -/,  sin  2fl>  =  2  sin  d>  cos  O, 
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eos24>  =  /-2  sin"  <£, 


tanf  /2 )  =  sin  /  ( /  +  cos  ) „ 

and  others  not  involving  inverses  of  sines  or  tangents  (i.c.s 
cotangents  or  cosecants)  because  Eq.  (Afi)  is  not  invertible. 

However,  the  exponential  matrix  lacks  one  important 
property — it  is  not  diagonalizable  via  a  similarity  transfor¬ 
mation,  the  usual  method  for  evaluating  the  exponential  of  a 
matrix.  Heren  the  exponential  is  evaluated  from  its  series  ex¬ 
pansion  (A4),  which  has  only  two  terms.  Moreover,  the  ar¬ 
gument  of  the  exponential  is  nor  dimensionless,  so  that  it 
must  be  understood  that  the  exponential  aces  on  column  vec¬ 
tors  of  the  form  {F.F')7,  where  the  lower  half  is  the  spatial 
derivative  of  the  upper  half  With  these  provisos,  the  use  of 
Eq.  (A4)  In  the  development  following  Eq.  (13)  is  internally 
consistent  and  achieves  notfltional  simplicity  and  transpar¬ 
ency. 


APPENDIX  B:  TRANSFER  MATRIX  METHOD 


The  EFA  boundary  conditions  require  the  continuity  of 


/  / 
l//]/ 2 


across  all  interfaces  (e.g.,  S.AfA  =  S fJ B) ,  where 


(BL) 


(B2) 


unless  (here  as  a  delta  function  discontinuity  in  the  Hamil¬ 
tonian  at  an  interlace,  e.g.,  ffXra^x)^t  in  which  case  an 
additional  phase  factor  comes  in.  For  example,  in  crossing 
the  BA  interface,  the  continuity  equation  reads 


W s=  **P( il ^  BA  )SJa  »  (B3) 

where 


exp  ( M's  J- 


I  0 

iHxyar}B  / 


Also,  at  AB  IF,  SjJA  =  e\p(-irt\4B)Sfif£  - 

The  second-order  Sehrbdinger  equation  is  changed  into 
the  first-order  equation 


df{i) 

m(z)/Kz)=-^? 

where  O  is  the  2NX2N  matrix 

0  /  \ 
-HTUMo-E) 


{B4jf 


(B5) 


where  K  is  diagonal,  Kif=  ,  i—  which  defines  a 

set  of  exponents  {A},  then 

exp(  /ft  j  )  =  P  txp{lKz )  P -■ 1 ,  (B8) 

where 

f  exptrATzJlj  -  exp{  lkp)£tj  .  (B9) 

integrating  Eq.  (B4)  across  one  period  from  —  b  to 
+  2a,  one  finds 

/1<{2fl+A)=exp(i2n>ja}/j,(A)^exp(/2fli<a} 

X£“l  e\p(-j(|>ffi,)SriCKp(f2nf&)/e(-6) 

=  exp(  i2il  ja)$j 1  exp(  -  «tp( 

X[5J!  (BIO) 

where  fAr#  >  SAiB  stand  for  quantities,  evaluated  at  the  At  B 
side  of  a  particular  interface;  in  Eq.  (BIO),  two  interfaces 
have  been  crossed.  Next,  employing  the  Bloch’s  theorem 

fA(2a+b)  =  tx${iqd)fA{-b)  (Bll) 

leads  the  secular  equation 

{[SA  txp(i2£lAa)SA  1  jexp(  ~i^BA ) 

^Kp{i2ilBb)Sg  1  ]exp(  -  expf 

Xfj-b)  =  0.  (B12) 

Define  matrix  M  [same  as  Eq.  (12)  in  the  text]  as 
S  =(SP)txpUKz)(SP)  ~] 

=  M^xpUKz)M~\  (B13) 

where 

M*SP  (B 14) 

Next,  Eq.  (BI4)  can  be  used  to  define  the  final  exponential 
matrix 

exp(  i  A  z )  —  M  cxp(  i  Kz )  M~ 1 ,  (B 1 5) 

so  that  the  secular  equation  [6 12)  becomes 

[  exp(  i  2  Aa  a)exp(  -  i  $  B A  )  txp<  i  2  A  sb )  exp(  i  3^  e) 

-eap(f?d)l/1y(-A>'0,  (B16) 

in  which  the  unknowns  are  energy  and  vector  fA{~b). 
Equation  (Bl6)  can  be  shown  to  be  the  same  as  the  Eq.  (16) 
in  the  text  and  can  be  used  to  derive  the  Kronig-Penney 
form,  Eq.  (30).  Clearly,  the  transfer  matrix  derivation  can  be 
generalized  to  any  number  of  layers  per  superlahice  period. 

APPENDIX  C;  ANALYTIC  MODELS 


In  regions  of  constant  composition,  O  is  constant,  so  that 
she  formal  solution  of  Eq.  (B4)  is^,6: 

/( z)  =  exp(  iffc  z }/{ 0 )  =  T(z)f(  0 )  «  (B  6) 

w‘here  /(0)  is  a  constant  of  integration  and  7"(z)^exp0'iiz) 
is  the  transfer  matrix.  If  P  diagonalizes  ft, 

F~'ttF=K3  (07) 


L  IF  eaupliag  of  HH  and  LH  bands  at  F 

For  the  coupling  of  the  HH  and  LII  bands,  the  4  X  4  prob¬ 
lem  at  r  decouples  into  two  identical  2X2  EFA  problems, 
so  that  all  hands  remain  doubly  degenerate.  Consider  then 
the  2x2  EFA  case  for  she  coupling  of  otherwise  uncoupled 
parabolic  HH  and  LH  bands  in  Eq.  (29).  For  parabolic 
bands12 
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tan  A  jfl  = 


/ 


0 
G 

0 


0 

0 

0 

fi2kLA  J2mL 


2/1  fifftff  /i 0 

U  2ALmL/h2kL 

0  0 

o  o 


(Cl) 


where  k'^_L=  ^2m'J{L{E-  V$)f&\  and 

is  the  valence  band  offset.  A  similar  equation  for  tan  A Bb  has 
kftfL=  For  symmetric  inter¬ 

faces,  the  interface  matrix  is  given  by 


exp(r^)- 


M 

o 

0 

\t 


0 

I 

-T 

0 


where 


T=  — (  — — — — j , 

i/5  \  2/n^j 


°\ 

0 

0 

1 


{CD 


(C3) 


bEy  degenerate  at  the  center  of  the  Brillouin  zone.  This  is 
consistent  with  the  symmetry  analysis  in  Sec.  SIT  D  for  sym¬ 
metric  $L%s, 


Z,  Anisotropic  IIII-SO  coupling  at  the  center  of  the  BZ 

Similarly  for  the  coupling  of  die  HH  and  SO  bands,  |f>  |) 
and  |]T  —  ^),  *>ne  finds 

(ffH  +  JtH  cot  kHb){  KS- ks  tan  ksb)  -  2{2T/h2  )\mHms )  =  0 

(CS) 

and 

( ks cotksb  +  KgXkff  tan klfb -  Kti)-\~2{2T/h2)2(mHms)=0, 

(C9) 


Altogether,  Eq.  (29)  for  a  symmetric  SL  at  I’  ean  be  solved 
to  yield  the  following  product  of  twro  determinants^ 

[A  BHmfyk*\klALlmAL+k{BL  tm* 

+  i2m1)2AHm^kii-{2m1)HAHmpkAfl)1=0 

m 

and 

[A  Lni'l /kf  +  B Lm^/kf  ]lkftA liftS u 

+  {2m2)ULmpkAil-{2T/ht)i{AtmtfkAL}1=0 

(C5) 

that  are  independent  of  the  sign  of  7'  and  are  real,  an  impor¬ 
tant  check  on  numerical  results. 

For  simplicity,  also  let  m^  L  be  continuous  across  inter¬ 
faces.  Then,  in  the  limit  of  a  quantum  well  (a -*■*), 
kHLA  tf.t-*  -  Xn,i  >  where  khl  =  /2ffiw(  V^Eyih1. 
Then,  Eq.  (C4)  becomes 

{kh+ kH cot  kHb){KL - kL  tan  kLb)  -  (2T/fi =  0. 

m 

Here,  the  first  and  the  second  factors  are  recognized  as  the 
conditions  for  odd  HH  and  even  LEI  bound  stales,  respec¬ 
tively  in  the  limit  of  no  LH  coupling.  Therefore,  the  cou¬ 
pling  connects  states  such  as  HH2  and  LHL. 

The  second  solution  yields 

(kj  cot kLb  +  KiK^F/  tan  kHb  -  kh)  +  (2 Tfh2)2{mHmL)  =  0, 

(C7) 

which  couples  odd  LI  I  and  even  HH  states,  e.g.,  Hill  and 
LH2-  This  equation  is  identical  to  that  derived  by  Ivchenko, 
Kaminskii,  and  Alciner.ifr 

The  equations  coupling  (he  other  two  states  |f,  —  j)  and 
\)  are  identical  to  the  ones  above,  so  that  all  bands  are  dou- 


where  =  'j2ms(E—  A)ft 2  and 

=  ^2ffTStA+  V0-E}fh\  The  factor  of  two  shows  that,  all 
other  things  being  equal,  the  IIII-SO  coupling  is  twice  as 
strong  as  the  HH-LH  coupling- 


3.  Anisotropic  HH-LH-SO  coupling  for  a  QW 

A  better  indication  of  the  combined  effect  of  die  HH- 
LH-SO  coupling  is  found  by  solving  a  model  of  a  QW  with 
the  coupling  of  ||a|>,  and  ||,“|).  The  Hamiltonian  at 

the  z  =  b  interface  is  given  by 


h1  d2 

2  rttjf  dz2 

iTS(z-b) 

V5T  S(z-b) 

—  iTS(z—b) 

ftJ  d2 

a. i2  d2 

■-il 

£ 

r-i 

! 

~2MdP 

ViTS{z—b) 

-ra. 

ih 2  d 2 

a1  d2 

~2 Tfd? 

2mg  dz2 

(CIO) 


where,  in  terms  of  Luttingcr  parameters  y, 

“2y2,  m%fmL=Y\  +  lyi>  m0/ffl5=yh  and  m^M 
—  2v^ y2 -  For  symmetric  interfaces,  at  j  =  T— ^  —  T. 

In  the  limit  of  large  spin-orbit  energy,  the  envelope  func- 
lion  for  the  symmetric  state  is  given  in  the  well  by 


fF„\  I  /hCOS  kfp  \ 

Fl  =  \/l  sin kL2-(imLiM)fs sin ksz  , 
\FS!  \  fssmksz  I 


(CM) 


and  in  the  right-hand  barrier  by 
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iFA 

Ft 

Fat 


I  (fHcoskffb)c\p[-  KH(z-b)]  \ 
ifi  sin  *L*)esp[  -  kl{z -ft)]  ~{imL !M) 
X  (/i  sin  k J*)exp[  -  ks{z- M] 

\  (/.TsinA-j6)cxp£- / 


by  the  LH-SO  mixing.  The  present  calculation  was  tested 
against  the  predictions  of  this  Appendix  using  mode]  QW's 
and  SL's. 

A  PPEN  DIX  D :  PHYSICAL  PARA M ETERS 


Z>b,  (Cl  2) 

where  the  continuity  of  the  envelope  function  at  interfaces 
has  been  imposed.  Here, 

A h,i  =  v'2 rfi //_ i Efh 2 *  A ^  —  i \'2 jtiA/i1?-  , 

kh,l=  v  -  to  f/,i(  /  A  ^ j 

+  A)/fi3,  and  the  "modified'1  spin-orbit  mass 
is  1  ip  =  1  /m  5  /  A/2 . 

integrating  the  Hamiltonian  equation  (CIO)  across  an  in¬ 
terface  yields  the  secular  equation 

h1 

t—  ( kH  tan  kHb  —  kh) 

2tnH 

T2  27*(1  +MtiftmL)2 

"" ^  _  ■ 

- — (itCoiAifi  +  tfi)  —  (KS+kscoiksb) 

(Cl  3) 

Using  the  physical  parameters  of  GaSb  from  Table  IX+  a 
comparison  of  the  HH-SO  model,  Eq.  (C9)  with  the  HH- 
LH-SO  model,  Eq.  (Cl 3),  shows  that  the  effect  of  the  SO 
hand  on  the  HH  levels  is  enhanced  by  the  factor 

( 1 -s-  MfV2mfy  Vji/ffrj=7_45  (Cl  4) 


APPENDIX  E:  LEVEL  SHUTS  FOR  VARIOUS 
COUPLINGS 

The  position  of  Ihe  HHl  level  at  T  was  examined 
using  four  modified  EFA  models,  Table  X.  Clearly,  the  shift 
of  the  HH1  band  docs  not  depend  linearly  on  the  strength 
of  the  Interface  interaction1  since  the  shift  of  the  HHl  band 
in  the  LHS  model  is  not  the  sum  of  the  shifts  in  the  HS  and 
HL  models  even  when  weighted  by  the  weights  of  the 
corresponding  bands.  In  fact,  the  total  shift  is  11.4  meV, 
while  the  sum  of  unweighted  shifts  is  2,6+ 1.6=  4.2 
meV  and  of  the  weighted  shifts  is  17.6  meVr  Observe  that  in 
a  QW,  there  is  a  great  density  of  hands  above  the  top  of  the 
QW  well,  which  must  be  taken  into  account  in  perturbative 
treatments  that  rely  on  the  completeness  of  unperturbed 
solutions. 

APPENDIX  F:  BOUND  STATES  OF  THE  LF  POTENTIAL 

Consider  a  hypothetical  type- 1 1  alignment  of  two  materi¬ 
als  A  and  B  with  otherwise  identical  band  parameters  but 
whose  interface,  nevertheless,  gives  rise  to  a  localized  inter¬ 
face  potential  coupling  parabolic  HH  and  LH  bands; 

H — // EFA  +HxraD3(z)®y  (FI) 

where  Hr fa  is  the  4X4  £c  ■  p  Hamiltonian  for  the  heavy  and 
light  holes  and  HXy  Ivchenko-Kamlnsky-Rossler  inter¬ 
face  coupling. 

The  EFA  equation  coupling  [I,  |)  and  \j?  -  manifolds 
becomes 


TABLE  IX.  Physical  parameters  used  as  input  m  Lhe  calculation 
(Ref.  63), 


Parameter  Symbol  GaSb  InAs 


^  d1 
~dP 

hh  i 


ft1/  / 

i/Z 


h1  tfi 


Conduction  band 

Ec  (meV) 

Overlap +80.3 

0 

ValeEtce  band 

E?  (meV) 

Overlap 

-410 

Band  gap 

Eg  (meV) 

813.3 

416 

Momentum  matrix 

Pz  (eV) 

25.  84 

22.04 

Spin-orbit  energy 

A  (meV) 

752 

3  SO 

Electron  mass 

mC 

0-04 1  mq 

0.023 

Luttinger 

7\ 

HE 

19.67 

parameter* 

72 

4.03 

8.37 

73 

5.26 

9.20 

Lattice  constant 

(A) 

6v0954 

6.0584 

Elastic  constants 

CL| 

9-00  X  10* 3 

8.33X10 

Cl2 

4.13X1011 

4.53X  10 

G44 

4.45  x  103] 

3-96  x  10 

Deformation 

6.85 

5.08 

potentials  (eV) 

n.79 

1.00 

by 

2.00 

-L8 

dv 

-4.8 

3-6 

\  2m0aD^^  2mLdP  / 

(F2) 

TABLE  X.  The  position  Of  the  HHl  level  ui  I’  calculated 
using  lour  modified  EEA  models  with  two  symmetric  IF’s  wtdi 
Wjfv=525  mcV.  In  ihe  2x2  H  model,  the  Hll  band  is  uncoupled; 
in  the  4X4  HS  model,  the  HH  and  SO  bands  are  coupled  by 
the  IFhs;  in  the  4X4  HL  model,  the  Ilf  I  and  LH  bands  are  coupled 
by  the  fF’s;  and  in  the  6  x  6  HLS  model,  the  interface  coupling*  arc 
HH-SO  rand  HH-LH,  so  that  LI  1-50  interface  coupling  is 
indirect. 


HHHmcV)  HH(%)  IH{%)  S0(%) 


2X2  M 

-50.78 

300 

0 

0 

4X4  HS 

-48.20 

99.86 

0 

0.14 

4X4  HL 

-49.17 

99.79 

0.21 

0 

6X6  HLS 

-39.40 

9B.57 

0,84 

0-59 
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where  the  energy  axis  is  up  (effective  masses  are  positive). 
Bound  solutions  are  sought  in  the  form 

FN=ffJ^xpi-KH\z\) 

?l  -fi  e*p{  -  kl  |z|  )*  (F3) 

where  khl  —  'J -2m  fJi jE/h 2  and  the  energy  is  negative.  The 
eigenvalue  equation  has  one  doubly  degenerate  hound  state 
at 


which  has  one  hound  solution  that  in  the  limit  of  large  SO 
splitting,  is  given  by 

rHjtmsqi  h2i  \* 

9fi*A  \2ma&J 

[£=-(0.008  meV^  for  GaSb]  (F6) 


E_  -  f  —  1 2 

[or  £  =  —  {2.97  meV)f2  with  GaSb  parameters]. 
Also,  for  a  coupled  H-  SO  system  (|j,  |)  and  |j,  -|))( 


(F4) 


ft2  <fi 

2m h  dP 


h2t 

2m0  oj 


4 


k2t 


\  2«io<Jo  "  3 


[2  a2  cP 


=£ 


so 


(FS) 


and,  in  the  limit  of  a  large  IF  potential  by 


r  ,  (  h2l.  |2*'WtfWSU 

\2mtaJ  3  ft3 


[£=-(2.58  me V)  t1  for  GaSb] .  ( F7) 


Therefore,  the  IF  potential  by  itself  can  bind  a  hole.  For  this 
reason,  the  problem  is  not  amenable  to  perturbative  treat¬ 
ments. 
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